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Abstract 

We discuss intensity effects in collisions between beams of optical photons from a high-power 



' laser and relativistic electrons. Our main focus are the modifications of the emission spectra due 

pi • to realistic finite-beam geometries. By carefully analysing the classical limit we precisely quantify 

, the distinction between strong-field QED Compton scattering and classical Thomson scattering. 

A purely classical, but fully covariant, calculation of the bremsstrahlung emitted by an electron in 
■ a plane wave laser field yields radiation into harmonics, as expected. This result is generalised to 

\^ , pulses of finite duration and explains the appearance of line broadening and harmonic substructure 

as an interference phenomenon. The ensuing numerical treatment confirms that strong focussing 
of the laser leads to a broad continuum while higher harmonics become visible only at moderate 
focussing, hence lower intensity. We present a scaling law for the backscattered photon spectral 
density which facilitates averaging over electron beam phase space. Finally, we propose a set 



X 

' of realistic parameters such that the observation of intensity induced spectral red-shift, higher 

harmonics, and their substructure, becomes feasible. 

PACS numbers: 12.20.Ds,41.60.-m, 

Keywords: high-intensity lasers, Compton scattering, Thomson scattering 
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I. INTRODUCTION 



Intensity (or strong-field) effects on QED scattering processes have been investigated since 



the 1960s following the invention of the laser. The pioneering studies considerec 



field pair creation [l| and the crossed process, electron photon scattering 2], [3, ^, ^,0,0. 



3oth strone 



where the use of laser beams has already been suggested. Since t 



len there has been a wealth 
for an 



0, Q, Q, Q, Q, Q 



of theoretical papers and we refer the reader to the reviews 
overview of the literature relevant for the present subject. From an experimental point of 
view the situation is less straightforward. There have only been a few clear-cut observations 
of intensity dependent effects. Probably the best known experiment is SLAC E-144 probing 
strong-field QED using a (by now moderately) intense laser beam in conjunction with high- 



energy electrons 



15|. Colliding a laser (L) of intensity 10^^ W/cm^ with the 46.6 GeV 



electron beam the observation of the nonlinear Compton scattering process 

e + £7l ^ e' + 7 (1) 

has been reported in [l^. Note that without an external field (here provided by the laser) 
an electron could not spontaneously emit photons as this is forbidden by energy momentum 
conservation. However, absorption of d. laser photons 7^^ induces the production of a high- 
energy (30 GeV) 7 quantum which thus takes away a large fraction of the incoming electron 
energy. This high-energy photon has then been used to produce electron-positron pairs 



171] via collision with the laser, employing the multi-photon Breit- Wheeler reaction 18| . 



7 + d'lL —>■ e'^ + e~ . Hence, using a high-energy setting with a large linac, SLAC E-144 has 
produced "matter from light" for the first time 

It is convenient to describe electron energy and laser intensity in terms of dimensionless 
parameters. The former is of course measured in terms of the relativistic gamma factor, 
7 = Ep/m = cosh(^ where ( denotes rapidity. A convenient measure of laser intensity is the 
dimensionless laser amplitude, 

eEX eE eE 
ao = = = , (2) 

with E being the root-mean-square electric field and X = 1/uj = 1/ mv the laser wavelength. 
Oq is thus a purely classical quantity, the energy gain of an electron across a wavelength 
measured in units of its rest mass. We mention in passing that quantum field theory imposes 
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an upper limit on uq [19| which becomes manifest upon writing 

eEl El 

ao = — - = — • (3) 

u Es V 

Here Es = m?/e = 1.3 x 10^*^ V/m is the QED critical field first discussed by Sauter [2Q] 
beyond which any laser becomes unstable due to pair creation from the vacuum via the 
Schwinger mechanism 2l|]. For an optical laser this implies the bound oq < l/i^ — 10®. 

The definition ([21) can be made explicitly Lorentz and gauge invariant 22] . When ao be- 
comes of order unity the quiver motion of the electron in the laser beam becomes relativistic. 
SLAG 13-144, for instance, had — 0.6 and 7 = 10^, i.e. low intensity and high energy for 
the purposes of this paper. 

From a classical (nonlinear optics) point of view the process ([1]) corresponds to the gener- 
ation of the ith harmonic in the 7 radiation spectrum. The production of higher harmonics 
has been observed in several experiments colliding laser and electron beams: low inte nsity 
laser photons (ao = 0.01) with low-energy ( ~ 1 keV)^ electrons from an electron gun 23| . 
ao = 2 photons with plasma electrons from a gas jet [24] and, more recently, sub-terawatt 
laser photons (ao = 0.35) with 60 MeV electrons from a linac 25|. Using linearly polarised 
photons the latter two papers have analysed the characteristic azimuthal intensity distri- 
butions confirming quadrupole and sextupole patterns for the second and third harmonics, 
respectively. The energy spectrum of the scattered radiation, to the best of our knowledge, 



has been measured only once, in an all-optical setup using laser accelerated electrons [26] . 
While this "table-top" setup is certainly attractive as it does not require a linac, the electron 
beam has a rather broad and random energy distribution which in turn is inherited by the 
scattered photons. As a result, the 7 spectrum recorded in 2^ is rather difficult to analyse 
theoretically. 

In this paper we discuss the prospects of experimentally analysing nonlinear Compton or 
Thomson scattering at a comparatively low centre-of-mass energy of the order of the electron 
mass, ^/s ^0.5 MeV but rather high laser intensity, / ^ 100 TW. Such an experiment is 
currently possible at the Forschungszentrum Dresden- Rossendorf (FZD), Germany, with the 
40 MeV linac ELBE ^ and the 100 TW laser DRACO [asj, so 7 = 80 and a maximum 
value of ao = 20. The linac provides a well-defined electron beam with high brilliance and 
low emittance so that a detailed study of intensity effects on the radiation spectrum should 
become feasible. In addition, new technology developed for ELBE enables the delivery of 
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"bunches" containing 1-10 electrons only. This could provide further new insights into the 
interaction of electrons with high intensity lasers. 

Our paper is organised as follows. To make the presentation self-contained Section [III 
briefly recapitulates the QED analysis of nonlinear Compton scattering. In Section IIIII this 
is compared with its classical limit (Thomson scattering). Section [IV] briefly discusses how 
to obtain the spectrum from a classical radiation point of view. Effects caused by the finite 
temporal and spatial extent of the laser and electron pulses are discussed in Section [Vl 
Section [Vlj contains our conclusions and summary. Some calculational details are relegated 
to the Appendices A and B. 



II. NONLINEAR COMPTON SCATTERING 



A. Kinematics 



Nonlinear Compton scattering, i.e. the processes ([T]) with ^ > 1, is treated in the text 
29 1 following [g] and has been reanalysed in 30 1. In this section we briefly collect the main 



results of the latter reference for later comparison, at the same time streamlining notation 
and normalisation conventions. Because of the large intensities involved one has to account 
for the possibility that many harmonic numbers ^ contribute to the total cross section. In 
other words, nonlinear Compton scattering is nonperturbative in both the electromagnetic 
coupling a and laser intensity Oq in that one has to sum over all individual processes ([1]). In 
terms of Feynman diagrams this is depicted in Fig. [H 

On the right-hand side we have displayed one term of the sum, the ^-photon contribution, 
a tree level diagram with dotted lines representing laser photons, solid lines electrons, and 
wiggly lines the emitted photons. The nonperturbative sum is depicted on the left-hand 
side, the double lines denoting effective or dressed electrons corresponding to the Volkov 



solution of the Dirac equation in a plane wave 



3l|. We are thus adopting a Furry picture 



32] where the interaction with the classical laser background is shuffled into the "free part" 
of the Hamiltonian having the Volkov states as its 'unperturbed' eigenstates. Schematically, 
the associated Volkov wave functions may be written as 

^p = e^^{iS[e]A]}T[e]A]up, (4) 

where Up is a free Dirac spinor, r[e;y4] a field dependent combination of Dirac matrices 




n7L 



+ 




FIG. 1: Feynman diagram for nonlinear Compton scattering (left-hand side) in terms of dressed 
electron lines (subscript V). This diagram may be expanded in an infinite series of standard QED 
diagrams in which the l-th term corresponds to the absorption of i laser photons jl (right-hand 
side). 



and 5'[e; A the classical Hamilton- Jacobi action functional for a charge e in a plane wave 
field A 291]. This suggests that the Volkov wavefunction is a WKB type solution of the 
Dirac equation. As expected from the Furry picture argument it contains all powers of the 
electromagnetic coupling e. 

Throughout this section we assume the laser photons to be circularly polarised which, 
for a harmonic plane wave, corresponds to the electromagnetic 4-potential^ 



A^^ = A'^{k ■ x) = a(e^ cos k ■ x + e2smk ■ x) 



(5) 



where fc^ 



= e,- ■ A; and e,- ■ Co 



'6ij. Assuming a head-on collision and summing over 



24 



25l | . On a more 



polarisation and spin states we thus expect the emission probabilities to be axially symmetric 
as there is no preferred direction other than the beam axis. Note that this is different for 
linear polarisation implying the azimuthal intensity patterns recorded in 
speculative level preferred directions along with azimuthal dependence may be induced by 
Lorentz violating physics such as noncommutative space-times js^. 

Suppressing electron spins the left-hand side of Fig. [1] may be analytically expressed as 



^ Linear polarisation is obtained by simply setting one of the e's to zero. The dimensionless intensities are 
then Uq = e^cP' jrv? and Oq = e^a} jlrv? for circular and linear polarisation, respectively. 



5 



the S-matrix element 

d'x^pix) (6) 

where e'^ is the polarisation vector of the emitted photon, k = {u, k) and k' = {u', k') are the 
momenta of laser and emitted photons, respectively, and p is the momentum of the incoming 
electrons before they actually enter the laser beam (assuming long but finite wave trains or 



adiabatic switching at infinite past and future 

As already mentioned in the introduction the analogous diagram in QED (with "naked" 
Dirac electrons) vanishes as one cannot satisfy p + k = k' with all three particles being 
on-shell. This is different in the presence of an external field which can transfer additional 
4-momentum. This statement seems obvious from the right-hand side of Fig. [H There is, 
however, an additional subtlety associated with the momentum assignment: in a plane wave 
field electrons acquire a quasi 4-momentum q reflecting the relativistic quiver motion. This 
may be seen by calculating the proper-time averaged 'Volkov current', jy = {^p'y^'^p) = 
q'^/p^ as in 29|] or by writing the Hamilton- Jacobi action in (jlj) as ^[e; A] = q-x + 0{{k-x)'^). 



In either case the resulting quasi momentum is the intensity dependent quantity 

g^p+^k, (7) 

2k ■ p 

and analogously for the outgoing electron momentum q' (with p replaced by p'). Working 
out the S-matrix element then yields a momentum conserving delta function with its 
support defined through 

P = q + ek = q' + k'. (8) 

The explicit dependence of quasi momenta on feeds through to the analogue of the 
Compton formula for the scattered frequency uj'. This is compactly expressed in terms of 
dimensionless scalar products upon writing the relevant four-momenta as 

k = Lon , k' = uj'n' , p = mu = 7m(l, f3) . (9) 

Note that only the velocity u transforms as a four-vector proper as only m is a world 
scalar. It is useful to measure frequencies in units of the electron mass, uj = um, uj' = u'm. 
Employing momentum conservation and the definitions (111) the (dimensionless) scattered 
frequency becomes 

= — 7 — — ; ■ ^ ^ 

n' ■ u + \ iu + — ^ ]n-n 
2n ■ u 



For the time being we will focus on head-on collisions where u = 7(1, —(3n) such that one 
can eliminate n' ■ u = n ■ u — jS'-f n ■ n' and ffTOl) simplifies to 
, Ivn ■ u 

^ ■ M + K£(ao) ^ ■ ri' 1 + K£(ao) (1 + cos6') 
Here, e"'' = 1/n ■ m = 7(1 — /5) is the usual Doppler factor, Q the scattering angle and 

K£(ao) = tv — sinh C + -Qq = — n ■ P/m (12) 

is the projection of the total momentum P = q+£k onto the optical axis, n = k/a;, measured 
in units of m. The vanishing of the latter, n - P = —vtiKu = 0, defines an intensity dependent 
centre-of-mass frame in which the scattered frequencies are precisely the harmonic multiples. 



z/' = iu 



For K£ < one is in the inverse Compton regime where the electrons transfer 
energy to the emitted photons, u'^ > iu, thus causing an overall blue-shift. This Doppler 
upshift, of course, is the physical basis for Compton generated X-rays. The maximum 
scattered frequency, i.e. the Compton edge, occurs upon backscattering {9 = 0) and is given 
by 

'^^.-x - i^2iue< + al • ^^^^ 
The presence of Oq in ( |T3i) leads to a red-shift of the n = 1 Compton edge compared to linear 
Compton scattering (oq 0) which, for 7^1, translates into the inequality 



Hence, if one is primarily interested in up-shifting the laser frequency (say, for X-ray gener- 
ation), the intensity should certainly not exceed unity. For large Oq the Doppler upshift 
may even be completely compensated due to the "stiffness" of the laser beam whereupon 
one leaves the inverse scattering regime. As shown in [30i], for 7^1, this happens when qq 
exceeds a critical value of 27. 

B. Emission rates and cross sections 

The S*- matrix element represented by the Feynman diagram of Fig. [T] (left-hand side) 



is readily translated into an emission rate 29|, l30l]. The differential rate (per volume and 



time) for emitting a single photon of frequency 00' = mv' in the process ([T]) is given by 

dWi e^m^ Ue 2 Ze{zi{x)) 



dx 16n gO "° (1 + x)^ 



(15) 
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where Ue is the density of incoming electrons. We have chosen a normahsation somewhat 
different from [s^ such that dWe now has dimensions of (i.e. particles per unit time and 
volume). The nontrivial part of the rate is encoded in the function 2^, 

Zi{z,) = -4 Jlize) +(2 + \jl,{ze) + JlM) - 2Jf{z,) \ . (16) 

The are the usual Bessel functions of the first kind depending on the invariant argument 



z,{x)=2^^^^J- (1--) (17) 



which is composed of two further invariants, namely 



^ = T-^^ yi = ^yi = r^, {0<x<ye). (18) 

ft ■ p 



Note that ze = when x acquires its minimal or maximal value. If we express the laser 
intensity in terms of (laser) photon density, n^, 

2 e^ni 

the differential rate f|T5l) may be written in a more symmetrical way, 

dx fcOgO (l + x)2 ' ^^^^ 

Tg = a/m ~ 3fm being the classical electron radius. Expressed in this form the rate is 



readily transformed into a cross section upon dividing by the (symmetric) flux factor 35| . 



where the last identity holds for a head-on collision. We thus end up with the differential 
cross section 

dx ' k-p{l + xy' ^ ^ 

which indeed has the correct dimensions of an area. Expanding Ze{z) for small one 
recovers the Klein-Nishina cross section jol. |29|^. 



^ Note that [2^ has shghtly different conventions: the densities ul and Ue are set to unity, and Gaussi, 
(rather than Heaviside-Lorentz) units are used which amounts to a reshuffling of factors of An. 
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For what follows we will need rates and cross sections in the lab frame where we assume 
a head-on collision between laser and electron beams. The scattered frequency is then given 
by (ITT]) and the invariants (ITHl) become 

!/'(! + cos _ 2ve^ 

^- eC-zy'(l + cos^) ' ^'"1+^- ^ 



This yields the following differential cross sections 

dap r^TT m? 



duj' n ■ P k ■ p 



Zi{zi) , (24) 
Zt{zt) . (25) 



The individual harmonic cross sections fl2^ are plotted in Fig. [2] (left panel). One clearly 
sees that the contribution of each harmonic has its own frequency range given by iuj < uj'^ < 
max' "^f- dlS]), with the individual supports overlapping to some extent. The contributions of 
higher harmonics, i > 1, are getting more and more suppressed in amplitude. These features 



are sufficient to guarantee convergence of the cross section summed over all harmonics 30|, 



da = Y^ dai . (26) 



=1 



The result of the summation (up to n = 500) is shown in Fig. [2] (right panel) for intensities 



Oo = 5 and Oq = 20. The figures basically coincide with those in 30|], but now we are able 
to state the absolute normalisation on the vertical axis. 

As already pointed out in [s^, two main features can be seen to arise. First, in line with 
our discussion at the end of the previous subsection, the linear Compton 'edge' is red-shifted 
by a factor of Qq ~ 400 (oq ~ 25) from about A'j'^u ~ 25 keV down to 0.06 keV (1 keV), i.e. 
from the hard to the soft X-ray regime when oq = 20. This is a rather drastic effect and it 



should be straightforward to verify experimentally [30|. Second, higher harmonics show up 
as additional peaks in the summed cross section with the peak heights decreasing rapidly 
with i. 

We conclude this section by pointing out that the observation of high harmonics may 
become feasible, albeit implicitly, using the total emission rate. For sufficiently large values 
of ao, very high harmonics, say even with i > 1000 have considerable integrated strength as 
shown in Fig. [3l For instance, when = 100 and i = 2000, the emission probability is 
only about two orders of magnitude below the value for the first harmonic, £ = 1. As Fig. [3] 
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FIG. 2: Left panel: Partial differential cross sections for the first few harmonics in nonhnear 
Compton back-scattering (head-on collision) for intensity oq = 20 (solid curve). Right panel: 
Harmonics summed up to £ = 500 for intensities ao = 20 (solid curve) and ao = 5 (dash-dotted 
line). The spectrum for linear Compton scattering (Klein-Nishina formula) is displayed in both 
panels for comparison (dotted curve). 

shows, summing such a large number of individual harmonics yields a rather smooth energy 
spectrum in total. 
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FIG. 3: The integrated emission probability = J dWg from (j20p as a function of the harmonic 
number i for various intensity parameters, ao = 1,5, 10, 20 and 100. 



For a more detailed analysis of nonlinear Compton scattering, in particular its lab frame 
signatures, we refer the reader to [30]. A complete discussion of polarisation effects has been 
given recently in jssl . 
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III. NONLINEAR THOMSON SCATTERING 



In this section we compare the results for nonhnear Compton scattering from the previous 



Section with t 
As stated in 



le results for nonlinear Thomson scattering considered, for instance, in |36l. l37l| . 



and further analysed in 30], the Thomson limit is obtained when the 



yeiao) = -^—^ < 1 . (27) 



invariant ye defined in (1181) becomes small. For head-on kinematics this is the statement 

2i z/e^ 
l + al 

Physically, the quantity ye represents the maximum (normalised) recoil of the electron during 
the scattering process. Hence, the Thomson limit amounts to neglecting the momentum 
transfer from the laser field to the electron. From fl27j) this is feasible if either ^ 27 and 
£ are sufficiently small or for large intensity, Oq, i.e. large photon density. The latter is, of 
course, consistent with approaching the classical limit. 

A closely related quantity measuring quantum nonlinear effects has also been introduced 



m 



61], namely 



ey/F^ k-p 1 
X= § — = — ^flo = 77?/i ao • (28) 

This has frequently being used since then, but as the last identity in fl25]) shows it is basically 
equivalent to (1271) . Following the example of 6|] we prefer to work with ye in what follows^. 
Finally, it is worth noting that ye is related to the usual Mandelstam variable 

s = iq + ikf = iq' + k'f = ml{l + ye) , (29) 

representing the total centre-of-mass energy squared. Hence, in the Thomson limit, one 
neglects the photon contributions to this energy, so that is the dominant energy scale, 
s ~ m\. In this sense, the Thomson limit is a low-energy limit. In contrast, the SLAG 
13-144 experiment [l^ (7 ^ 10^, oq — 0.6, hence yi ~ 1) has probed the genuine Compton 
(or quantum) regime. 

In what follows we want to find explicit relations between the general Compton expres- 
sions and their classical (Thomson) limit. To this end we try to separate off the quantum 
corrections from the purely classical results. We begin with the momentum projection (IT^ 
which may be rewritten as 

Ke = kq + ^?/£(1 + ao)e~^ , (30) 



The second equality in (pS)) follows from the gauge invariant definition of ao, see e.g. [22 1 
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where kq is obtained by setting i = 0. Replacing ne — ^ /to in the scattered frequency (fTTj) 
for head-on colhsions straightforwardly yields the (nonlinear) Thomson limit, 

iun ■ u £u „ , 



'^'^'^^ n-u + Kon-n' 1 + Koe-'^{l + cos 6*) ^'^^ 



(31) 



This suggests that the general formula for arbitrary collision geometry is obtained by setting 
£ = in the denominator of (fTOj) such that the scattered frequencies 

^^;,Th = = -2 (32) 

n- u -\ n ■ n 

2n ■ u 

are indeed integer multiples of a fundamental frequency u'^. At this point it is instructive 
to compare with the low intensity ("linear") limit (or Thomson limit proper) where i = 1, 
flo — > and (l32l) condenses to 

u' = z/ . (33) 

This is the Doppler shift in disguise upon noting that for a head-on collision and backscat- 

tering n ■ u = 7(1 + /5) = and n' ■ u = 7(1 — [3) = e~^. 

Expressin g th e invariant x from ( fTSj) in terms of the scattering angle 9 it becomes explicitly 

^-dependent [301], x = = ixi, with 

^ z/e-^(l + cosg) 
' l + /«oe-qi + cos^) ■ ^ ^ 

Comparing with ( |3Ti) one finds the relation 

= • (35) 

As xe is bounded by the Thomson limit implies X£ — > so that Compton and Thomson 
expressions should generically differ by terms of order as in ( |35l) . This is consistent with 
the findings in jscj] and confirmed by the numerical comparison of Fig. HI left panel. 

Let us move on to the emission rates. Again, we try to isolate all terms dependent on xe- 
Following [3^ we define the ratio r = Xi/yi with < r < 1 and rewrite Z£ from (fT7|) as 



z, = izi = 2e ^ Vr(l - r) , (0 < ;2i < 1) . (36) 



'I + a2 



It is important to take the Thomson limits, Xi —>■ and yi —>■ 0, in such a way that the ratio 
r stays fixed as a result of which ze remains unchanged. We may therefore decompose the 
spectral function (fT6!) into a classical (Thomson) part and an xi dependent correction, 



^2 



Z(iz() = M^e) + —^2e{ze) , (37) 

i -p X£ 
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where (suppressing the overall argument z^) 

£^ = J^_^ + J|^i - 2 J| , and = -U^/al + 2£.e . (38) 
This yields the Thomson limit of the cross section ( |25l) . 

where only the classical part, of contributes. Factoring ^£ out from fl25|) we may 
explicitly calculate the Compton to Thomson ratio, 

Again, Compton and Thomson results differ by terms of order xi. A numerical comparison 
is displayed in Fig. IH right panel. 




FIG. 4: Left panel: Comparison of classical (Thomson) and QED (Compton) results for the 
scattered frequency in the backscattering direction, = 0, with bv'^ = xh ~ ^'i plotted as a 
function of y^. Solid line: xi/{l + xi), cf. (|35|) : dotted line: Thomson limit (|31|) . Right Panel: 
Normalised integrated partial cross section ai as a function of the electron recoil ye normalised to 
Thomson limit, cf. (j40p . Choosing i = 1 the curves for oq — > and ao = 5 up to 20 are on top of 
each other. 

For our set of parameters one finds ?/£ ~ 2 x 10~^i which is very small unless one considers 
extremely high harmonics {i > 10^). Thus, according to (11CT|) . Compton and Thomson results 
differ by approximately 10~^ for low harmonics. This suggests, for these parameters, that a 
purely classical calculation will yield a very good approximation for the Compton scattering 
results. This will be further corroborated in the following Section. 
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IV. CLASSICAL RADIATION 



Having considered the classical limit of Compton scattering it is worthwhile comparing 
with a purely classical calculation. This will later be used to investigate various finite size 
effects on the spectrum. We are particularly interested in a fully covariant description of 
the radiation emitted as bremsstrahlung by a charge moving in a laser field such as given 
by ([5]). It is most convenient (and, as an additional bonus, also maintains explicit gauge 
invariance) to start from the field strength tensor [2^ 



F^'' = F,{k ■ x) fr , fr = ^<-rf4. ^ = 1,2, (41) 

with = k'^/uj as before and ef the two transverse polarisation vectors introduced for the 
circularly polarised gauge field which corresponds to the choice'' 

Fi{k ■ x) = —au sin(A; ■ x) , F2{k ■ x) = aw cos{k ■ x) , (42) 

such that F = auo represents the magnitude of electric and magnetic field. 



The Lorentz equation of motion for a charge in a plane wave has first 



been solved by Taub in 1948 



and quite elegantly) 



(for other early treatments see |36|, 



39, 



40|). A modern 



covariant analysis has recently appeared [2^ on which the following remarks are based. For 
the following it is useful to absorb charge e and mass m into the field variables i.e. to replace 
{e/m)F^^ — i> F^^^ and (e/m)A^ A'^ such that the rescaled A^^ is dimensionless^ . The 
Lorentz equation of motion then reads 

p'^ = F^,(x(r))j9% (43) 

and is, in general, a nonlinear differential equation as the field strength, F'^'^, depends on 
the trajectory x(r) to be solved for. Thus, normally, one cannot expect to find an analytic 
solution. For a plane wave, however, where F'^" = F^"{k ■ x) there is a sufficient number 
of conserved quantities such that the system becomes integrable. Most important is the 
constancy of longitudinal momentum in time, k-p = mVt = const, which defines a frequency 
fl = uin ■ u such that k ■ x can be traded for proper time via 

k ■ X = Qt = ujn-UT = ujn-UoT, (44) 



^ To obtain linear polarisation just set one of the F's to zero. 

^ The intensity parameter hence becomes Oq = —{A'^) = a} for circular and a§ — a? /2 for linear polarisation, 
respectively. 

14 



adopting the initial conditions k ■ x(0) = and m(0) = uq. Integrability is now obvious 
as F'^'^ = F'^'^{Qt) which turns fH3|) into a hnear equation. All amplitude functions be- 
come functions solely of proper time. Integration of is now straightforwardly done via 
exponentiation, u = exp(J dr F\un. The null field f l41l) is nilpotent of order two, i.e. all 



powers higher than two vanish 22]. Thus the exponential series gets truncated after the 



term quadratic in field strength leaving us with 

OA . 7,„ _ /|2 

MA'(r) = - + 2 , = A^(t) . (45) 

Note that we have assumed the initial condition A'^(O) = 0, otherwise one has to replace 
^ — ^4^(0). It is reassuring to directly check the conservation law n-u = ii-uq which 

holds as n ■ A = (Lorenz or light-cone gauge) and = 0. 
Another r integration of ( H5l) then yields the orbit 

a;^(r) = x^(0) + r< - / dr' A^" + f dr' {2A ■ Uo - A^) n'' . (46) 

Jo 2n • uq Jq 

The first two terms are obviously initial conditions, followed by a transverse part along A'^ 
and a longitudinal contribution proportional to = k^/uj. 



The calculation of the c 
trodynamics (see e.g. [41 



assical radiation spectrum may be found in most texts on elec- 



42]) though not necessarily in a covariant manner. One has to 
determine the 4-momentum of the radiation field using Poynting's theorem and the re- 
tarded potentials, Ajtt = □"^j^]ret- A particularly compact covariant expression is given by 



the Fourier integral 43] 



P' = -l -^,sgn{k')6{e)k^j{k).j{-k), (47) 
where j^{k) is the Fourier transform of the electron current, 

j/^(x) = e [ dX^S^x - X(r)) = e [ dTu''{T)5\x - X(r)) , u" = dX^/dr . (48) 



Accordingly, (1471) expresses the radiation 4-momentum in terms of the electron trajectory, 
X = X{t), encoded in the current ( HHI) . The zero component of ( H7I) . of course, yields the 
radiated energy which, after performing the k^ integration, we write as 

dJdVl J ' J , k' = n') , (49) 
duj'dil 
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where we have reinstated primes for scattered momentum components. The integrand is 
the spectral density describing the number of photons radiated per unit frequency per unit 
sohd angle, 

^2 A/- / 

Thus, to determine the spectral density all we need to know is the Fourier transform j^(fc') of 
the current ( HHI) which depends both on orbit position and velocity, and u'^, respectively. 
Employing the continuity equation, k' ■ j{k') = 0, one may eliminate j° such that (|5U]) turns 
into 

p(^',n') = ^|n'xj(fc')l'>0. (51) 

In summary, a determination of the classical radiation spectrum in an external field amounts 
to solving the Lorentz force equation in this field for the trajectory, and hence, to find the 
current j^(x) of fj48l) . Its Fourier transform, compactly written as 

jA^(fc') = ej rfrM^(r)e-^^'-^(") , (52) 

then yields the spectrum via the radiation formula (l50l) . The Fourier integral ( l52l) may be 



(partly) evaluated using a trick of Schwinger's [4^, assuming periodic r dependence of the 



gauge field as in ([5]) and fH2l) . In this case the four- velocity is periodic as well, 

M^(r + e'T) = n^(r) , T = 27r/n , f integer , (53) 
but for the orbit one finds 

x^(r + e'T) = e'Tw^" + x^(r) , (54) 

where w'^ is the velocity averaged over one period of proper time and hence proportional to 
the quasi-momentum i^, = mw^. Note that the dependence on £' has been separated off 
in the first term of ( |5^ . For (almost) periodic functions /(r) it makes sense to decompose 
the integral over all r into a sum of integrals over all periods, i.e. 

I dTf{T)= J2 [ drf{T + e'T). (55) 

Applying this to the current (152!) and, once again, separating off the C' dependent pieces one 
finds the expression 



ik') = e ^ exp |27rzf ^ ■ uq + ■ n'^ | j^{uj', n') , (56) 
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where we have defined the integral over a single period T, 



J dTU^'^r) exp i^-iu' n' ■ Xq - iuj' j rfr'n' ■ n(r')| , (57) 



which is independent of i'. Thus, in f l56|) all dependence on i' has again been factored off such 
that the sum over all periods can be evaluated using Poisson resummation, exp{27Tii' h) = 
Yle^i^ ~ ^ result, we obtain a "delta comb" for the current, 

nk') =uj[J2 ^(^' - M) Jt(M, n') , (58) 
e>o 

where the multiples of uj[ = mv[, derived from the exponent in fl571) . define the harmonic 
frequencies v[ = C.v[ which precisely coincide with the Thomson limit (132|) upon identifying 
Mo with the asymptotic velocity u = p/m of the scattering process^. Mod-squaring our 



answer we thus conclude with Schwinger [4^ that periodic motion (induced by periodic fields 
of infinite spatio-temporal extension) leads to a line spectrum of radiation into harmonics 
labelled by index i. 

What can we expect to happen when we restrict the fields to have finite temporal duration, 
i.e. for pulses? In general, an analytic treatment will be difficult, but there is one particular 
case which is reasonably straightforward and nevertheless yields the basic physics involved. 



This is the case of a rectangular pulse (see also [45||), where we just cut off the r integral in 
fl52|) such that —tq/2 < r < tq/2. Accordingly, the sum over £' in fl56|) only extends from 
— to A^, where we assume that our pulse contains 2N periods, tq/T = 2N. The finite sum 
can nevertheless be evaluated with the result that fl58|) gets replaced by 

uj' 

sm{2N + 1)tt— 

j%{k') = jt^iuj', n') . (59) 

sm TT — 

The prefactor composed of the ratio of two sines is nothing but the diffraction pattern 
obtained when light passes 2N + 1 slits. It has the usual maxima at u' /uj[ = i integer but 
also A^ — 1 additional peaks of lower amplitude between adjacent integers. Thus, in our 
language, each harmonic is accompanied by a substructure of A^ — 1 = tq/2T — 1 secondary 
peaks. Furthermore, there is of course line broadening, as the (formal) zero-width limit of 
the delta comb is only achieved for infinite A^, hence infinite temporal extent. 



^ For head-on collisions we recover (131 
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These findings should qualitatively also hold for more realistic (smooth) pulse shapes 
such as Gaussians or power laws. We will explore this in more detail in the next section. 

Before we come to that we conclude the general discussion with a remark on the relevance 
of radiation damping. In principle, the Lorentz equation fl43l) is only valid approximately 
as it neglects the back-reaction of the radiation field on the electron orbit (see 46(] for an 
illuminating recent discussion). To incorporate the latter the Lorentz equation is superseded 



by the Landau-Lifshitz one 
lasers 



42l | which was recently utilised in the context of high-intensity 



48| . For our parameters, the radiation reaction parameter defined there becomes 

4 



R 



10" 



(60) 



We have checked that for our electron beam and laser pulse specifications the longitudinal 
motion is barely altered. In particular, radiation damping modifies the change in energy after 
collision by just 0.58 %. Accordingly, the backscattered photon spectrum will change at the 
sub-percent level at most, and we can safely neglect back-reaction effects. For more detailed 



discussions of radiation damping in a laser context the reader is referred to [Sg, |47|, |48|, |49 |. 



V. BEAM SHAPE EFFECTS 



It has become clear that treating the laser beam as a plane wave of infinite spatial and 
temporal extent according to (I5l) and ( l42l) is an idealisation, the validity of which has to 
be checked. An infinite plane wave should be a fair approximation if the electrons only 
probe the central region of the laser focus and if the reaction time is small compared to 
pulse duration a. To assess the feasibility of these assumptions and the size of possible 
modifications we will study finite size effects by first considering pulsed plane waves and, in 
a second step, allowing for a transverse intensity profile in addition. 



A. Finite temporal duration 

A pulsed plane wave is obtained upon multiplying the field strength F^^" from fHTl) cor- 
responding to the infinite wave ([5]) with an envelope factor of width a such that (H2l) gets 
replaced by 

Fi{k ■ x) = —ga{k ■ x) F sin(/c ■ x) , F2{k ■ x) = g„{k ■ x) F cos(/c ■ x) , (61) 
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with some suitable envelope factor and F = auo. Note that this does not spoil the fact 
that the associated gauge potential (in radiation gauge) still represents a plane wave 
solution of Maxwell's equations, as k"^ = uj'^n^ = implies □A^(/c ■ x) = 0. By writing the 
pulse width as o" = Qtq with a pulse duration Tq in proper time, we have 

g^{k ■ x) = gir/To) . (62) 



Orbit velocity and trajectory are still given by (H5l) and fH6l) if the pulse factors are included 
in the definition of the gauge potential A^{t). Hence, a useful analytic expression for the 
orbit can be obtained whenever the r integrals over A and A'^ can be evaluated which, of 
cou.e. cucally depend, on *e pulse shape fu.c.o I. has bee. suggested use 
the "solitonic" pulse 

5,(r) = sech(r/ro) , (63) 

which, being the generating function of Euler numbers, does not lend itself to straightforward 
integrations. A particularly simple case, however, is obtained by considering crossed fields 
with F'^'^ constant, hence A^ linear in r, 

F^"" = Ff^'', A" = aQr . (64) 

Superimposing a power-law pulse shape, the amplitudes become time-dependent, e.g. 

F oQt 

such that the orbit coefficient functions appearing in become 

^ dr'A^(r') = aQry (^^I + t^t^ -1^ , (66) 
r dr' A\t') = -a2^]V3 f- - arctan . (67) 

Jo V^O ToJ 

As a result we note that within the pulse (r ^ tq) the coefficients (!66l) and (!67|) are quadratic 
and cubic in r, respectively, hence behave as for infinite crossed fields. Outside the pulse, 
i.e. when it has passed by (r ^ tq), the orbits become modified and approach free motion 
linear in r. 

As we have seen in the previous section, for any particular choice of amplitudes one has to 
insert the velocity (Il5|l and the orbit (H6|) into (!52|) and evaluate the radiation formula (150|) 
to analyse the infiuence of the finite temporal extent encoded in the pulse shape g{T). For 
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FIG. 5: Normalised spectral density p = p/tq as a function of normalised frequency to' je^^io for 
a = 0.7 (left panel) and a = 1.4 (right panel). Upper (lower) panels correspond to a laser pulse 
duration of Tq = 25 (Tq = 100 fs). The full vertical (red) lines denote the nonlinear Compton 
edges for the ideal case of an infinite plane wave laser field, while dotted (grey) lines represent 
the respective linear Compton edges. For pulsed fields the spectral density covers the whole range 
between the two Compton edges. 

a rectangular pulse we saw secondary peaks appearing in a way reminiscent of a diffraction 
pattern. 

For smooth pulse shapes g the situation is slightly different. Adopting backscattering 
kinematics for simplicity and the circularly polarised pulse (ET]1 . the electron current 
becomes 



-e J drA^r) exp < —iu' | n' ■ uqT 




n ■ uq 



(68) 



where the dependence on proper time r in the exponent is strongly nonlinear due to non- 
vanishing a ~ 1 and the presence of the smooth pulse envelope g. 



As a numerical example, we have chosen the pulse shape (jHS]) advocated for in 50| and 
determined the radiation spectrum for laser amplitudes a = 0.7 and a = 1.4 and pulse 
durations^ Tq = 25 fs and 100 fs assuming circular polarisation as before. In Fig. [5] we have 
plotted the normalised spectral densities p = p/To as a function of the normalised frequency 
uo' /e^'^uo. Looking at the spectra one observes Nj. sub-peaks within the first harmonic signal 



The pulse duration Tq in the lah frame is defined by setting J d{k ■ x) ga{k ■ x) = wTq max {g^). For the 
solitonic pulse Tq is related to the proper time pulse duration tq via Tq = ire^TQ . 
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resulting from the nonlinear r-dependent modulation in the exponent of (l68i) . Radiation 
generated at different times r, and therefore at different effective laser intensities a?g'^{T) 
with effective Compton edges u' /e^'^oj = (1 + a^5'^(r))~\ interferes thus generating the 
pattern of sub-peaks seen in Fig. O There are important differences between smooth pulses 
and the interference for box-shaped fiat top pulses which were discussed in f l59p and below. 
Here, the interference pattern is due to the nonlinearity of the process (a ~ 1) together 
with the non-trivial envelope function g of the laser and it disappears in the limit a ^ 0. 
The smooth pulse sub-peaks are a sign of chirp in the emitted X-ray pulse. This however, 
is not the case for box-shaped pulses, where the side-peaks just reflect the discontinuous 
pulse shape and resemble the interference pattern of a single slit. These differences will be 
further investigated in Appendix [B] where the time structure of the scattered radiation is 
discussed. Summarising we conclude that sub-peaks in the harmonic signals are a purely 
nonlinear effect. 

It has been known for a while that temporal modulations strongly affect the spectral 
densities, with the main effect being the additional oscillatory substructures 45|, |47|, l51 |. 
Our analysis above gives a rather simple explanation of this phenomenon. Interestingly, 
a similar pattern (with a similar explanation) has been observed recently for the rates 



associated with laser induced pair creation 52| which is obtained from Thomson/Compton 
scattering via crossing. 

We have found that the number N-r of subsidiary peaks within the first harmonic scales 
linearly with the pulse duration Tq and the intensity o? according to the empirical formula 

AT, = 0.24ro[fs]a^ (69) 

Reversing the arguments leading to (1^ suggests the interesting possibility to actually de- 
termine the intensity of a laser pulse by counting the number of sub-peaks within the first 
harmonic if the pulse duration is known. 

To the best of our knowledge, the harmonic substructure has never been observed exper- 
imentally. This may be due to the fact that most certainly the sub-peaks will get smeared 
out by a variety of mechanisms. These include effects such as (i) the ponderomotive force in 
focused beams, (ii) the influence of transverse intensity profiles as well as (iii) contributions 
due to the phase space distribution of the electron beam, in particular its energy spread 
A7/7 and transverse beam emittance e. These effects will be studied in detail in the next 
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subsections. In particular, we will show how these effects can be minimised to allow for 
a possible experimental verification of the individual subpeaks with a 100 TW class laser 
system. 



B. Spatial intensity profile 



Let us briefly turn to the effects associated with a transverse intensity profile of the laser 
beam. This is taken into account by choosing the complex vector potential 



A = ae'^ib, z) g^jik ■ x) e 



ik-x 



(70) 



with a (linear or circular) polarisation vector^ 



e e 



{e^, e.y} for linear polarisation, 

{e+, e„} = {e^; + iey, — iey] for circular polarisation. 



(71) 



and the transverse distance h = [x^ + y^)^^^. The paraxial approximation of the wave 
equation then yields the following equation for the transverse profile ^ {s^, 



1 -i 



9a 
9cj 



(72) 



with the transverse Laplacian Aj^ = 9^ + 9^ and primes denoting derivatives with respect 



io k ■ X. The standard Gaussian beam solution (with = 1) reads 531] 



^o(p,2) 



Wo 



exp 



, , , . ^n I exp<2arctan — 
w{z) { w{zy J \ zr 2R{z 

with focal spot radius ("waist") wq, Rayleigh length zr = WqUj/2, beam radius w{z) = 
wo{l + z^ / ZrY^"^ and the curvature of the wave fronts R{z) = z{l + z\/ z^). 

In order for the Gaussian profile \E'o to be an approximate solution of the paraxial equa- 
tion fl72l) . the temporal profile function g„{k ■ x) = g^r/To) has to satisfy \g/g\ ^ ^ [50 1 
The Gaussian gi^r) = exp{—T^ /2tq) does not have this property, since \g/g\ ~ r is un- 



ci; 6^ 



(73) 



bounded. For a hyperbolic secan t, q (r 



\g/g\ = tanh(r/ro)/ro is bounded jSOl. l54|. 



sech (r/rg), the condition can be fulfilled, as 



Note that the use of unnormahsed polarisation vectors {e+,e_} is required here to be consistent with 
the normahsation of in (ISl) . 
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For a strongly focused beam such that A = wq/zr = 0(1), the intensity profile \&o will 
have to be corrected, since the fields derived from \[^o solve Maxwell's equations up to terms 
of O(A^) 54 . |55| . Although these corrections to \E'o a^re crucially important in some cases 



19|, 



561], we did not find them relevant for the electron trajectories or the emitted photon 



spectrum associated with the (almost) head-on collisions studied in this paper. 

An important effect of the transverse beam profile is the ponderomotive force Fp = 
— VmA^/2 pushing the electrons away from regions of high intensity as they gain trans- 
verse momentum Qj. Accordingly, electrons with /3 || before the scattering will leave 
the interaction region under an angle ctout with respect to the z-axis. For fixed total 
pulse energy Wtot and pulse length Tq, the magnitude of the ponderomotive force scales 
as Fp oc Wtot/ (^o^o) oc max (oout)- This means that the ponderomotive force leads to mea- 
surable effects only for very small waist size wq. To quantify the effect we simulate electron 
trajectories corresponding to head-on collisions assuming a laser pulse of 3 J energy and 
To = 20 fs for different impact parameters b (cf. Fig. [6]). For wq = 5 yum the maximum 
deflection angle for 40 MeV electrons is about 1 mrad and for 10 /im it is roughly one order 
of magnitude lower. This is small compared to the typical angular scale of the emitted 
radiation which is of order I/7 ~ 12 mrad. 




20 30 
b [/iin] 

FIG. 6: Effect of the ponderomotive force deflecting an electron with impact parameter b, initially 
in perfect alignment with the laser pulse (zero injection angle), after collision ejected with an angle 
Oout with respect to the beam axis. Laser parameters: energy 3 J, pulse duration To = 20 fs, 
different focal radii wq. 



Let us move on to our numerical simulations, choosing a 3 J laser pulse described by 
either a circularly or hnearly polarised Gaussian beam of pulse duration Tq = 20 fs and 
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wavelength 800 nm colliding head-on with a dilute electron beam. The spectral density for 
Ne electrons is calculated as an incoherent superposition of the individual emission rates pi 
according to 



Figure [7] shows the spectrum p{uj' , 6) in the plane = for a tightly focused laser with 
Wo = 5 /im and a corresponding Rayleigh length of zr = 100 pm. The peak values of the 
normalised amplitude, a = 8.66 and a = 6.12 for linear and circular polarisation, respec- 
tively, are clearly in the nonlinear regime, the associated nonlinear Compton edge being 
approximately 1 keV, cf. flTTl) . The electron bunch — in this subsection we use cold electron 
bunches — is modeled by a Gaussian with a bunch length of 1 ps and a transverse beam 
size of r;, = 5 pm. The spectral density, normalised to one electron, is shown as a function 
of the energy uj' for different scattering angles, 6 = 0,5 and 10 mrad. Due to the strong 
field gradients of the laser in both transverse and longitudinal directions, the spectrum is 
extremely broad for the parameter values chosen. The individual harmonics are not visi- 
ble as the individual spectral lines are overlapping. It seems fair to describe the scattered 
photons broad continuum. 

In comparison. Fig. [8] shows the backscattered spectrum for the same pulse (3 J, Tq = 20 
fs), but with a larger focal radius of Wq = 50 fim corresponding to a Rayleigh length of 
Zr = 10 mm. The maximum value of the normalised amplitude is a = 0.866 (a = 0.612) for 
linear (circular) polarisation, which corresponds to a nonlinear Compton edge of 27.5 keV. 
The electron beam has the same transverse beam size ri, and bunch length as before (Fig. [7j), 
so that the electron bunch exclusively probes the very centre of the focus where the laser 
intensity is almost constant. In fact, at = 5 /im = wq/10 the field intensity is only 
1% lower than at the centre of the focus. In contrast to Fig. [71 the harmonics are now 
well-separated and clearly visible. 

C. A Scaling Law for the Spectral Density 

For the following we assume that the laser waist size is much larger than the electron 
beam radius, Wq 3> r^, so that the electron beam interacts only with the centre of the laser 
focus. In this case the laser may be reasonably described by a pulsed plane wave in the 




(74) 
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FIG. 7: Normalised spectral density p{oj' ^ 0)/N(, of the scattered radiation as a function of frequency 
u)' and different scattering angles, = 0, 5, 10 mrad (from top to bottom). Assumes scenario: head- 
on collision of a dilute electron bunch {N(, = 10000) with a strongly focused {wq = 5 /um) laser 
pulse (3 J, 20 fs). Left panel: circular polarisation. Right panel: linear polarisation. 
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FIG. 8: Same as Fig. [7] but for wq = 50 /xm. 
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relevant interaction region 



We may thus concentrate our attention on energy spread and emittance.^ According to 



59( 1 the former may be modelled by a distribution of the Lorentz factor 7 of width A7, 



centred at 70, while the latter measures the transverse phase space volume of the beam via 
the correlator 

e. = loPo ^/WWy^Wf - loPo nAa. (75) 

The expectation values (• ■ ■ ) refer to the transverse phase space distribution of the electron 
ensemble, usually taken to be Gaussian as well. In addition, we have defined the normalised 
transverse momentum, ^ = \px/Pz\ = tana ^ a, which basically coincides with the injection 
angle a with respect to the beam axis (chosen as the z direction). Obviously, while our 
previous considerations were assuming a head-on collisions {(3 ■ n = —(3), we now have to 
consider sideways injection allowing for a small angle a relative to the laser beam axis (see 
Fig. [9]), i.e. /3 ■ n = —[3 cos a. In this case we have to use the general relation (|TTi) for the 




z 



FIG. 9: Sketch of beam geometry. The initial electron velocity is denoted by u, while n (n') is the 
direction of the incoming (outgoing) photon. 

scattered frequency depending on three scalar products. 



n ■ n 



1 + cos 6 , 



n ■ u = 7(1 + /? cos a) , 
n ■ u = 7(1 — [3 cos a) . 



(76) 
(77) 
(78) 



^ In fact, things may be turned around by using Thomson scattering as a diagnostic tool to measure electron 
beam parameters 58|. 
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Here, 9 is the scattering angle between k' and — k, a is the angle of incidence between 
laser and electron beams (p and — k) and a' the angle between p and k'. For a head- 
on collision, a = and a' = 9. Allowing for the possibility of linear polarisation we 
choose ei = so that laser direction and polarisation define the xz plane. This introduces 
another azimuthal angle, if, for the electron momentum, p = 7m/3. We thus have in general 
(5 = /?(sina cos ip, sin a sin ip, cos a), n' = (sin 9 cos0, sin 9 sin 0, cos 6*) and 

cos a' = cos(6' — a) — [cos(6' — a) — cos{9 + a)] sin^ ^ ■ C^^) 

As we now have a distinguished vector transverse to the beam, namely p_|_ = {px,Py), we are 
breaking axial symmetry and thus expect the radiation to develop an azimuthal dependence 
on (p. The same is known to happen for linear polarisation (even for head-on collisions, 



a = [37j) and has been used to detect higher harmonics [24 1. 

To proceed we need the dependence of the spectral density fl50l) on the initial conditions, 
in particular the initial 4- velocity mq. Clearly, this enters via the current (l52l) according to 



f{k') = e rfrM^(r;Mo)e-^'=' "("'"o) . (80) 



Writing uq = 7o(l,/3o) with /3q = /3o(sin ao cos y^o, sin ao sin (^o, cos ao) we see that the spec- 
tral density will depend on the initial electron energy 70 (via /5o) and the angles of incidence, 
ao and ipo, p = p{uj', n'; 70, Oq, ^o)- 

How will this change with initial conditions? To address this question we view such a 
change as being due to a Lorentz transformation, A : mq — ^ ma = Auq, where A is defined 
in Appendix |X1 Under this assumption one can derive a scaling formula that relates the 
spectral density for the transformed initial conditions to the original one. Choosing the 
initial kinematics of a head-on collision (oq = 0, ipo = 0) for simplicity and u\ = ma(7, a, ^) 
we find 

p{uj', n'; 7, a, if) = M(7, a, ip; 70, 0, 0) p{uj' / h, n ; 70, 0, 0) . (81) 
Here, we have defined a rescaling factor for the scattered frequency 

h ^ 4^ (82) 

with z/'(m) as in ( l32l) and a 'transition function' which depends only on kinematic quantities 

z/'(ua) {n-UoY J l-4s + 4s^^ linear polarisation 
u [Uo) [n ■ ua) I 1 _ 2s + circular polarisation 
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where we have introduced 

uj [n ■ ua) |e ■ ua\ 

Here, e = (0, e) denotes a single, real or complex, polarisation vector describing linear 
or circular polarisation fITT]) . Averaging M corresponding to linear polarisation over the 
azimuthal angle ip one naturally recovers M for circular polarisation. The first factor in 
( l83l) is related to the Jacobian of a Lorentz transformation for going from head-on to side- 
injection geometry and accounts for the fact that the radiation is peaked in the forward 
direction of the electron. The derivation of the scaling law (18T|) is deferred to Appendix |Al 
The main virtue of the scaling property fl8T]) obeyed by the spectral density p is its use 
in calculating what is called the 'warm spectral density' p^yi^u'). The latter is defined as the 
expectation value of the ('cold') spectral density p(7, a, (f ) taken in the initial ensemble of 
Ne electrons, characterised by the normalised distribution f{'-f,a,ip) of initial energies and 
injection angles. Thus we have 



= I d-fdadip /(7, a, ip) M{-f, a, cp; 70, 0, 0) p{uj'/h; 70, 0, 0) . (85) 



As pointed out in 59|], in this way one avoids to perform a summation over an ensemble of 
test-particles of the form fITil) . The distribution function f{'y,a,ip) in (185!) is taken to be a 
product of a Gaussian in 7, a x-distribution (with 2 degrees of freedom) in a and a uniform 
distribution in (p, 



a 



2 



/(7,«,^) = ^mUa) = (^^^yiT^exp |-^(^| (^exp ("^(^ j (86) 

as we assume the electron beam to be axial symmetric with both transverse components u^o 
and Uyo normally distributed and uncorrelated. In order to resolve the substructures in an 
experiment, one has to fine-tune several parameters. First of all it seems sensible to only 
have a small number of subpeaks (say less then five), so that they can be clearly resolved. 
This is achieved by adopting values of a = 0.7 and Tq = 25 fs (these values correspond to the 
upper left panel of Fig. [5]). It is important to have good control of both the energy spread 
and the emittance. The energy spread tends to smear out the small subpeaks, in particular 
at the high energy end of the spectrum, whereas emittance affects the whole range. 
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FIG. 10: The warm spectral density for ELBE 



27l | parameters (left panel), LWFA electrons 



Ml 



(centre panel) and for our proposed parameters (right panel, see Tab.|T]). For comparison the cold 
spectral density is shown also in each plot as grey dotted curve. The smearing of the harmonic 
sub-peaks due to finite emittance and energy spread is clearly visible. Very low values for both the 
energy and angular spread are needed for a possible observation of the sub-peaks. 

A common source for ultrarelativistic electron beams are linacs like the ELBE accelerator 
at the FZD [27]. It is capable of producing electron bunches with a very low energy spread 
of A7/70 = 10~^ and transverse emittance of about = 1.5 mm mrad. On the other hand, 
new laser based acceleration schemes like laser wake field acceleration (LWFA) report the 
production of electron bunches with A7/70 = 3.5% and angular divergence of Aa = 0.68 



mrad at 70 = 400 GOj. The numerical results for the warm spectral densities using both 



ELBE and LWFA electron beam parameters are shown in Fig. [10] which also includes a 
comparison with the respective cold spectral densities. Clearly, the harmonic sub-peaks 
are smeared out for ELBE parameters and even more so for LWFA. For ELBE it is the 
emittance which is too high for the observation of the sub-peaks. A LWFA electron bunch, 
on the other hand, has too large an energy spread — despite its low angular divergence. We 
conclude that in order to resolve the sub-peaks both energy spread and emittance need to 
be sufficiently small. In Tab. [T] we list a suitable set of parameters which allow to observe 
the harmonic sub-peaks with a 100 TW laser. For petawatt lasers, the strong constraint on 
the emittance may be relaxed because they are capable to achieve nonlinear peak intensities 
> 1 across larger spot sizes wq such that larger electron beam radii can be tolerated. 
The essential quantity in fISB]) is Aa oc e/vb and not the emittance itself. 
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TABLE I: Beam parameters required for observing the substructure in the fundamental harmonic 
of nonUnear Thomson scattering employing a 100 TW laser. 



parameter 



proposed value 



laser frequency 
laser amplitude 
laser pulse length 
laser focal radius 
electron beam radius 
electron bunch length 
electron energy spread 
electron transverse emittance 



a;[eV] 
a 

To [fs] 
Wo [/im] 

Lb 

A7/70 

£x [mmmrad] 



1.5 
0.7 

25 
50 

5 

<^ zn = 10 mm 
0.001 
0.7 



VI. DISCUSSION AND SUMMARY 



Backscattering of an optical laser beam by relativistic electrons has become important as a 
tunable source of X-rays (see e.g. 6l|)- At low laser intensities (oq ^ 1) the relevant physics 
is adequately described in terms of the conventional classical or QED treatment (Thomson 
or Compton scattering, respectively). At high intensities (qq ^ 1), however, one enters the 
relativistic nonlinear regime where multi-photon processes become important. These require 
a strong-field QED approach adopting a Furry picture where Volkov electrons dressed by 
the external field replace the ordinary QED electrons. As a result, the backscattered photon 
spectrum deviates significantly from the one corresponding to the standard Klein-Nishina 
formula for (linear) Compton scattering. By the correspondence principle, one expects that 
large photon numbers should allow for a classical description of the laser beam. Indeed, we 
find that classical Thomson scattering (i.e. bremsstrahlung by an accelerated charge in the 
external laser field) yields the same answer (e.g. for unpolarised cross sections) as strong- 
field QED at leading order in the small recoil parameter, Xi. For an optical 100 TW laser 
and 40 MeV electrons, such as in operation at Forschungszentrum Dresden- Rossendorf, one 
has Xi ~ 10~^ £, where i is the number of laser photons involved. The sharp decrease of 
emission probabilities with photon number i corroborates the validity of working in the 
Thomson limit. In particular, as we have shown, effects due to finite size and realistic shape 
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of both laser and electron beams may be easily addressed within this framework. 

Our main focus was the search for suitable conditions allowing to experimentally observe 
the signatures of the non-linear, multi-photon processes in question. In particular, we have 
found that the rich substructure in the spectral density of the first harmonic strongly de- 
pends on the combination of short pulses and high intensity. To quantitatively assess this 
dependence, we established a simple scaling law for the spectral density by means of which 
we could estimate the effects of electron emittance and energy spread on this substructure. 

Turning to experimental prospects we conclude that the ELBE/DRACO constellation at 
Forschungszentrum Dresden-Rossendorf offers the possibility to clearly detect the red-shift 



of the Compton edge (as already argued in 30||) as well as the generation of higher har- 



monics (i.e. their intensity distribution and angular dependence) as long as the laser pulse 
is not focussed too strongly (oq = 0(1)) and of sufficiently short duration. Higher inten- 
sities, accomplished by strong focussing, lead to a near-continuous backscattered radiation 
spectrum reaching far beyond 100 keV in the ultraviolet. This spectrum emerges from the 
superposition of many higher harmonics, modified by the temporal and spatial variations of 
the laser pulse. This renders a clear-cut verification of high-intensity signals rather difficult. 

In summary we have presented realistic results for the Thomson/ Compton backscattering 
spectra of optical laser photons by relativistic electrons assuming an experimental set-up that 
can be realised at facilities already in operation. The subtle interplay between beam and 
intensity parameters may be fine-tuned in such a way that the observation of intensity effects 
due to the increased effective electron mass seems feasible for the first time. 
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APPENDIX A: DERIVATION OF THE SCALING LAW 

The scaling law (IHT!) relates the spectral densities p for different geometries, in particular 
for different initial electron velocities, Uq. It is useful to describe such a change in geometry 
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as being due to a Lorentz transformation, A : Mq — -u^ = h.^^u" , where A is a Lorentz trans- 
formation composed of transverse rotations Kr and a boost A^ which may be parameterised 
by 

/ 1 \ 

COS a + sin^ — COS a) — sin cos ^9(1 — cosa) cos(y9sina 
— sin(y9 cosv9(l — cosa) cosa + cos^ ^^(l — cosa) siny^sina 



Ab(C) 



y — cos if sin a 

^ cosh C — sinh C, ^ 
10 
1 

Y — sinh COO cosh C, ) 



sm ip sm a 



cos a 



(Al) 



(A2) 



for a rotation with angle a around the axis = (— sin(^, cos(^, 0) perpendicular to the z- 
axis and a boost along the z-axis with rapidity C = cosh~^ 7, changing the electron's energy, 
respectively. We construct A as follows: First rotate uo(7o, ao, '/'o) parallel to the 2;-axis, 
then apply the boost changing the electron's energy from 70 to 7 and finally rotate to the 
new direction a, 

A = kda, (^)AB(7)A-i(7o)A-i(ao, <^o). (A3) 

According to ( l80l) the spectral density p depends on the initial conditions through the orbit 
x(r; Uo) and its velocity u{t;uq). In what follows, the direction of observation n' is kept 
fixed. Let us consider the simplest case first, namely backscattering and a head-on collision, 
implying a change only in the electron energy, 70 to 7 (as a = ip = 0, Afi = 1). In this case, 
the spectral density obeys the simple scaling relation 



p(w',7) = p{uj'/h,-fo) 



(A4) 



with a rescaled frequency u'/h where h = z/|('Ua)/z/((-Uo) and u'^ as in (!32|) . 

If we also allow for a change in the direction of uq, the height of the spectral peak will 
certainly change, since the radiated intensity is peaked in the direction (3 of the electron. 
Thus, we have to use the modified ansatz p{uj'; u\) = M{u\, Uo)p{uj'/h] uq) with a transition 
function M which we calculate in what follows, for low intensities (a^ <^ 1) and arbitrary 



initial value Uq, but strictly staying within the backscattering geometry, i.e. n' 



-n, 
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n' ■ n = 2. Linearising the orbit expressions (H5|) and (H6|) in the gauge field A, 



n ■ uq 



x''{t;uo) = 4 + / dr' A>'{t') + I dr'^^""'^ ' 



n ■ Uq 



(A5) 
(A6) 





the three- vector part of the electron current (l80l) becomes, to 0{A), 



}{uj') = e dr 



A ■ Uq Lu' I t I a / i\ 

n — A — 2i Uq / dr Air ) ■ Uq 

k-uo n-uo J 





(A7) 



It is convenient to adopt a complex-valued vector potential A'^ = (0, A) with A 
aee*^'^5'(r/ro)) with Qt = k ■ x and the polarisation vectors e as in (ITTIl . 
The inner integral in ( ]A7I) yields, after an integration by parts. 



-A'^{T){l + 0{l/nTo)) 



(A8) 



turn ■ Uq 

For sufficiently long pulses, Qtq ^ 1, the second term can be neglected. Within these 
approximations, the result for the electron current reads 



-ea J dr g{T)e~ 



■i{Lj' n' ■uo—Lon-uo)T 



Jo, 



(A9) 



where 



k ■ Uq 

Using flSTl) the spectral density becomes 



e ■ uo uje-uo 
}o = 7— n + e - 2-7— ^uq. 



uj{n ■ u^y 



\g{uj -uj'{u/iy'))\'\n' x jo 



(AlO) 



(All) 



IGvr^ (n ■ Mq)^ 

with z/' from fl33|) . Here and in the following g denotes the Fourier transform of the envelope 
function g, and 



|n' X joP = |ep + 4t2|e ■ UoP(u2 - (n • Uo)^) - 4t|e ■ Uo|^ 



(A12) 



with t = uo' /uj {n ■ Uq)"^. The form of g in (lAlip lends support to the scaling behaviour of 
the frequency adopted in flA4l) . On changing the geometry, uq — > ma, and replacing the 
frequency, uj' uj' = hu', with h = u' {ua) / u' {uq) , the function g remains invariant. 
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To determine the transition function M, we chose the special head-on geometry charac- 
terised by uq = 7o(l, 0, 0, /5o) as a reference, and a second, different geometry characterised 
by Ma = Amq or, equivalently, the injection energy 7 and the injection angles a and (p 
(cf. Fig. [H]), which yields 

p(a;';Mo) oj' (n • ma)^ |n' x jo(7o, 0, 0)|2 
Evaluating (lA12p with some explicit polarisation vectors one eventually arrives at (1831) and 
f l84p ■ To approximately include the non-linear case when a > 1, we finally substitute the 
nonlinear scattered frequency fl32|) for u' fl33l) in the definition of h, i.e. 

v'^v[ = . (A14) 

n' ■ Uq + n' ■ n 

2n ■ uo 

In the nonlinear regime, the scaling law perfectly describes changes in the electron's initial 
energy. However, changes in the angles of incidence, in particular a, are rendered accurate 
only for a ^ 1. Nevertheless, for the purposes of this paper, the scaling law is sufficient to 
account for the typical angular divergence of the electron beams. 

APPENDIX B: TIME STRUCTURE OF THE BACKSCATTERED PULSE 

To complete our analysis of effects due to finite pulse duration let us briefly comment 
on the temporal structure of the backscattered pulse. The Fourier transform F{t) of the 
electron current (IA7p . 

00 

F(t) = J dco' jico') e"^'' (Bl) 
—00 

provides information on the time structure of the scattered pulse in the lab frame. If the 
envelope function g in the current fl68|) is chosen as a 'solitonic' pulse as in fl63|) . the inner 
integral may be evaluated analytically with the result 

j dr'g^r') = rotanh|^-^^ . (B2) 

Thus, (IBip becomes 

F(t) = —e (irA(nr, t/to) 5 ( t — n' ■ uqt -\ — Xotanh — 

J \ ^ n-uo tq. 

= -eaexp(C) f(t) (B3) 
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FIG. 11: Time structure of the scattered X-ray pulse for a circularly polarised laser with 'solitonic' 
pulse shape (a; = 1 eV, a = 10 eV~^, = 10). Left panel: Electron current f defined in (IB4p as a 
function of time t. Red (light) line: low intensity (a ^ 0); black line: high intensity (a = 3). Right 
panel: Normalised time dependent frequency as a function of time t. Red (light) and black lines 
as before. Additional green curve: intermediate intensity (a = 1). Dotted lines: Thomson limit 
{i = 1, v ^ 0) of scattered frequencies p^ . 



with 



m 



A{w{t)t,t/a{t)) 



1 + 



(B4) 



where we defined a time dependent effective width cr(t) via the transcendental equation 

t 



t 2 , ^ 
+ a tanh ■ 



(B5) 



a{t) a{t) To ' 

and a time dependent frequency zxj{t) = ujc'^Tq/ a{t). By construction, the product of width 
and frequency is constant, w(t)a{t) = cue^TQ. The inversion of (]B5|) has to be done numer- 
ically. The results are shown in the left panel of Fig. [11] for two different values of a and 
for cij = 1 eV, To = leV"^, which corresponds to a pulse length of Tq ^ 20 fs in the lab 
frame with exp((^) = 10. The main features are (i) an increase of the scattered pulse length 
for larger a and (ii) a double chirp of the backscattered pulse due to the time dependent 
frequency zu(t) - the frequency decreases towards its minimum at the centre of the pulse 
and then increases. 

The chirp in the backscattered signal is a combined effect of the non-linear interaction and 
the non-trivial envelope function g with finite pulse length a. Both features are required for 
the chirp to be present. Hence, if either a ^ or a ^ oo the chirp vanishes. For box-shaped. 



35 



flat-top envelope functions, the time dependent frequency zu{t) wiU be discontinuous at the 
beginning of the pulse, jumping from its linear value e^^cu to its nonlinear value, e^^a;/(l-|-a^), 
and back again at the end of the pulse. In the intermediate regime zu{t) will stay constant 
so that there will be no chirp for this particular case. 
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